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Abstract 

The Tavis-Cummings model of N two-level atoms interacting with a single resonant mode is extended to various 
cases of off-resonance, initial photon densities, and atom number N. For stimulated absorption and small numbers of 
mean photon number a very unexpected results showing well defined rise and fall of oscillation is found. 


Background; 

The time development of the ensemble averages of E~ and E~E^ for photon number densities 
representing the eoherent and thermal states [i] was given in 1965 but only for the resonant ease and short 
times [2]. Eberly [3] provided the long term development for the ensemble average o^E~E^ showing the 
periodie eollapse and revival of oseillations in that average. An elegant empirieal demonstration of the 
quantum eollapse and revival was observed in 1987, when a single Rydberg atom in a single mode of an 
eleetromagnetie field in a supereondueting eavity was investigated [4]. Sinee that time, numerous papers 
have been published on the subjeet of quantum radiation in interaetion with atoms. More reeently, results 
for simulated emission in a single mode were given for both resonanee and non-resonanee for various 
initial photon density distributions [5]. Ini968, the exaet solutions to what has beeome known as the 
Tavis-Cummings Model were presented [6,7], extending to N atoms the basie model of one atom 
originally proposed in 1963 to study the relationship between the quantum theory of radiation and the 
semi-elassieal theory in deseribing spontaneous emission [8]. 

In the following, we will apply the solutions of the Tavis-Cummings model for N two-level 
moleeules (N-TLMs) as basis states to find expressions for the time development of the ensemble 
averages of f and E~E^. With the general solution in hand, speeifie eases of interest will be presented 
here in more detail for various photon distributions previously developed [5]. 

Ensemble Averages for E~ and E~E^ 

The ensemble average of the field operators E (t) and E E^(t) may be found in the 
usual way [9] 


{E (t)> = “ (“) ^(” + l)^/^<n|p/(t)|n + 1), 

^ n=0 


and 
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(E“E+(t)) = |-| ^n{n\pf(it)\n), (lb) 
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where only one mode of the field is exeited, y is the eomplex eoupling eonstant, and p 
is the dipole moment of the TLM with whieh the field is interaeting. The elements of 


the field density matrix are given by the traee over the TLM states 
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where 
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and p(t) is given by a unitary transformation of the density operator at time to =0 where 
it is assumed that the N-TLMs and radiation field are not interaeting! Therefore 

p(t) = U(t)p(0)U-Ht), (4) 

y(t) = (5) 

and where H is given [6,7]^ by 

H = (a^a + R^) + a^a —— - xaR^ — (b) 

Sinee the system is non-interaeting at time zero, the density operator is a direet produet 
of the field part and N-TLM part of the system. 

p(0) = Ptlm ® Pf. (7) 


Equation (2) ean be expanded direetly in terms of the orthonormal eigenveetors of 
the Hamiltonian above assuming that all the TLMs are at equivalent mode 
positions, namely the |r,e,j> states whieh are expressed in [^], 

c+r 

\r,c,i >= ^ > \r,c — n>. (8) 
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In expression (8), r and e are good (eonserved) quantum numbers of the Hamiltonian where ^ - and 

e=n+m with —r<m<r. Then 

{n\pf(t)\n') ^ ^P(r) ^ ^ {n\{r,m\e~^^^\r',c',j') 

r,m ^ ^ r',c',j'r",c",J" 

X {r',c',j'\p(0')\r",c",j"){r'',c'',j"\e'-^^\r,m)\n'). 

^ Reference 7 contains the full development for both resonance and non-resonance while reference 6 contains only 
the resonance results for which the photon energy o) is equal to the energy separation Q for the TLM. 





With this and a fair amount of effort [10], equations (la) and (lb) ean be written as^ 
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As is seen, these equations are quite eomplieated and to use the general expressions would be prohibitive 
exeept for the simplest eases. This eompletes the formulation of the general expressions. We only eonsider 
two speeial eases next. 


Special cases 

Equation (10) ean be simplified eonsiderably by eonsidering speeial eases for the initial TLM distribution. 
Towards that end, we eonsider eases for simulated emission and absorption where all the TLMs are 
initially in the up state or down state. 

All TLMs initially in the up state-Simulated Emission 

For this case ^ = ^ andP(r) = 1. Note that r is the cooperation number, n the photon number and 
c=n+r. Further 
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Equation (10) becomes 
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^ The factors and |^| seen in eq. (10) differ from those same constants in reference [2] by a 
factor of 2, i.e. y is 2y in that reference. We ignore the difference here. 
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In addition, the density matriees (c-|py^(0)|c-) and (c-|py^(0)|c + 1-) only have values for 


c — ^ > 0. By ehanging variables and using the orthogonality relationships for the following 
expressions ean be found [101: 
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where n is the mean number of photons and 
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Equation (14) is the final general solution for 5i(n,Af,7t) (no assumptions about resonanee or 
non-resonanee) and we have used the alternate expression for the effeetive 
eigenvalues of referenee [6], namely X=c-\K\q and multiplied by Q to obtain the 
eorreet units even though h is still set to unity. Note that n\K\t = yt. Using the 
notation introdueed by referenee [2], the ensemble average for the field {E~(t)) ean be 
written as 
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and we have taken advantage of the real nature of the ^4^. The relative tuning parameter is a measure of 
non-resonanee and is given by (3 = Simplifieation of ^2 is diffieult exeept for resonanee. We will 

noteonsider {E~(t)) or S 2 (n,N,yt) further in this paper. Speeifie examples for N=1 
through 4 for resonanee and N=1 and 2 for non-resonanee were provided in referenee 
[ 10 ]. 

Numerical evaluation of 5^ (n, N, yt) 

The exaet ealeulation for this expression is performed using eq. (14). In order to perform these 
ealeulations we make use of the photon number densities developed in referenee [5], namely Eqs. 17-28 
in Table 1. 


Table 1: Photon Number Distributions 
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^ Roy J. Glauber, Coherent and Incoherent States of the Radiation Field, Physical Review, Vol. 131, #6, Sept. 1963, 
p. 2766. 

^ M. S. Kim, F. A. M. de Oliveira, and P. L. Knight; Properties of squeezed number states and squeezed thermal 
states; Physical Review A, Vol. 40, # 5, Sept. 1, 1989 

^ J. J. Gong and P. K. Aravind, Expansion coefficients of a squeezed coherent state in the number state basis. The 
American Journal of Physics, Vol. 58, Issue 10 Oct. 1990 

































Table 1: Photon Number Distributions 
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See reference [5] for definitions of various functions used in Eqs. 17-28. In the following examples, we 
will not focus on numerical solutions for N=1 TLM since those cases were examined in some detail for 
the various photon densities in reference [5]. We will provide some numerical solutions for the photon 
densities represented by Eqs. 17-18 for N=10 as a reminder of what we considered in references [10, 
15 ]. The photon density representative by Eq. 19 will not be considered. All of the examples provided in 
the Table 2 below are for 10 TLMs since calculation for larger values takes considerable time. Each row 
of the table contains the photon density on the left and the value of Si(n,N,Yt) as a function of time 
yt on the right. Again n is the mean photon number and N the number of TLMs. In some cases we 
provide only resonant cases. For others we provide only non-resonant cases such as for the thermal 
state and squeezed thermal state. It is noted that numerous examples could be provided for various 
values of coherent photon number, thermal photon number, squeezing parameter, phase, and initial 
photon number and for various off resonant parameters. The range of possible examples is too large to 
explore in this paper, thus I provide the Mathematica notebook used to perform the calculations. 


^ A. Vourdas, Superposition of Squeezed Coherent States with Thermal Light, Phy. Rev A Vol. 34, #4, Oct. 1986, p. 
2366. Note that the expressions for the mean and variance were not given in this reference. In fact, the authors have 
not found this expression for this case in any reference. 

^ Paulina Marian and Tudor A. Marian, Squeezed States with Thermal Noise. I Photon-Number Statistics, Phy. Rev. 
A, Vol. 47, #5, May 1993, p. 4474. 

^ F. A. M. de Oliveira, M. S. Kim, P. L. Knight and V. Buzek, Properties of displaced number states, Phy. Rev. A, 
Vol. 41,#5,p2645 

^ P. Krai, Displaced and Squeezed Fock states. Journal of Modern Optics, Vol. 37, #5, p889, 1990. 


































I did explore a small range of coherent photon number for the coherent state represented by Eq. 17 
since this was an easily calculated example. When the number of photons is equal to or smaller than the 
number of TLMs, we see from the examples below that the photon number is insufficient to drive the 
TLMs to a state of equal up and down TLMs and that the reoccurring oscillations rapidly approach 
random oscillation. 

It is believed that the rise and fall of oscillations seen for the non-resonant cases is due to the distortion 
of bases states away from the more or less harmonic state forms exhibited for the resonant cases. This 
distortion can be extreme to the extent of favoring only a few of the at the lower or upper values 


of n. 


In Table 3 below, we provide examples which examine the effect of increasing photon number for the 
resonant coherent case represented by Eq. 17. Two sets are shown, 1 for 1TLM and the other for 10 
TLMs. As in Table 2, Table 3 is in a row format with the photon density on the left and the value of 
Si(n, N,yt} as a function of time yt on the right. The rows alternate with one row for 1 TLM and the 
second for 10 TLMs. We see that as the photon number increases greatly beyond the number of TLMs, 
that the rise and fall becomes considerably sharper and further apart. This is attributed to the linearity 
of the differences of eigenvalues in Eq. 14. 

By examining the exact solution found in Ref. [5] displayed in Eq. 29 for one TLM, we see that as n 
increases we can approximate Eq. 29 by the following equation 30 where A is the standard deviation in 
photon number. 


OO 



29 


n=0 



30 


As is easily seen, the square root can be expanded as a constant plus a term linear in the difference 
between n and n. 

Note that some of the examples shown within Table 3 are repeats of examples in Table 2 for ease of 
review. Further, for the last example in the table, time was only taken out to 800 rather than 2500 since 
a time estimate of between 4 days and 600 hours for computation was made. The second special case 
will be considered next. 

All TLMs initially in the down state-Simulated Absorption 

For this case, r = ^, m = — r, and P(r} = 1. Starting with eq. (10) and applying the various restrictions 
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This expression is very similar to eq. (14) but with some differenees sinee the upper limits on the sums 
over j, j’ and p depend on n and the eigenveetors y4^’^'''are summed downward from the highest term 
rather than upwards from the lowest term. It can also be seen that for N TLMs, that there are N-1 unique 
terms, one each for every n term less than N. It is only forn > A1 that all the terms look the same. 


Numerical evaluation of N, yt) 

The calculations for this case appear to be no more difficult that for all the TLMs up, although the 
precision for which the calculations are performed must be raised to prevent round off errors or division 
by zero. As with the TLMs all up, all the photon distributions provided in Table 1 can be used for 
example calculations of 54 (n, N, yt). Examples for all cases will not be provided but can be obtained 
using the supplied Mathematica notebook. Instead, the examples provided will at first focus on the use 
of the coherent distribution for various numbers of coherent photons and a comparison of the results 
for 54 (h, Al,yt) and Si(n,N,yt) (all TLMs up) for the same number of TLMs. Some interesting results 
will be seen. For a large number of coherent photons compared to the number of TLMs, the results will 
appear quite similar, although with differences. However, when the number of coherent photons 
approaches and then becomes smaller than the number of TLMs, striking differences between the 2 
functions will occur. The results are presented in Table 4 . There are nine sections in Table 4. Each 
section has both a left and right column and is 3 rows deep. The left and right columns are read 
vertically with the top row of each section the graphic for the photon density. The second row is the 
case for all TLMs up (5i(h, N, yt)). The bottom row is for all TLMs down (S 4 (n, N, yt)). Only resonant 
cases are provided although the reader may consider other cases using the supplied program. Section 1 
contains the 2 cases of 10 coherent photons ( Eq. 17 in Table 1) and 5 TLMs in the left column while the 
right contains the case of 100 coherent photons and STLMs. As mentioned above, the case of mean 
photon much larger than the TLM number results in very similar and 54 . In section 2, the cases of 
smaller mean photon number compared to the number of TLMs are presented. In fact the left column 
has a mean photon number of 10 vs. 100 TLMs while the right column contains the case of 25 mean 
coherent photons vs. 100 TLMs. The cases for the TLMs initially up is much more chaotic than the case 
of all TLMs down. Again the photon density was calculated from Eq. 17 in Table 1. The results in the left 
column for all TLMs down is very surprising in that a well-defined rise and fall in oscillations is seen. This 
behavior prompted a continued calculation in Section 3. Both left and right columns were for 100 TLMs. 
The left column was calculated for 1 coherent photon while the left was for 100 coherent photons and 
again the photon density is provided in Eq. 17, Table 1. It is again seen that there is a well-defined rise 
and fall of oscillation for the left column with all TLMs down. That is not the case for the right column 
but another unexpected results is seen. Namely that the average number of absorbed photons is more 
than have the number of mean photons. This is surprising since in all previous calculations for TLMs in 





the up state, the maximum number of mean emitted photons was never more than Vi the number of 
TLMs and often much smaller (the larger the number of mean photons, the closer to Vi the mean 
number of emitted photons for the resonant case). In the case here, the mean number of absorbed 
photons is approaching 60. (Admittedly, this number is now dependent on photon number rather than 
TLM number). We will discuss this somewhat later. 

Based on the previous results, it was decided to continue the investigation of a small mean number of 
photons vs. the number of TLMs but with the other photon densities provided in Table 1. Section 4 
contains the results for the Thermal photon density in Eq. 18 was used for both left and right columns. 
The mean photon number was 1 for the left column while the right had a mean photon number of 10. A 
TLM number of 100 was again used. Again we see a very surprising results, namely the expected chaotic 
behavior for the TLMs up but a well-defined oscillation for the TLMs down. This was not expected since 
all previous calculations with the Thermal Photon density always displayed chaotic behavior unless 
strong non-resonance was applied. 

Sections 5, 6, 7, 8 and 9 contain similar examples for the other photon densities in Table 1. For all the 
remaining cases the number of TLMs was constrained to 50 to limit the time it took to generate the 
figures (This time is listed in the figure call out). Section 5 has the results for the photon densities 
defined by Mixed Coherent and Thermal State (Eq. 20) in the left column and the Squeezed Coherent 
State (Eq. 24) in the right. The calculation parameters for the photon densities are shown in the call 
outs. Section 6 has the results for the photon densities defined by Squeezed Thermal State (Eq. ^) in 
the left column and Squeezed Vacuum State (Eq. 21) in the right. Section 7 has the results for the 
Squeezed Fock State (Eq. 22) in the left column and the DSTS (Eq. 2^ in the right. Section 8 has the 
results for the Mixed Squeezed Coherent plus Thermal Noise State (Eq. 25) in the left column and the 
Displaced Number State (Eq. 27) in the right. The final section 9, contains only the Squeezed Displaced 
Number State photon density (Eq. ^). In nearly all of these cases, the time development of is rather 
chaotic although if a few cases there is some coherence for a short time. However, again for all TLMs 
down we see rather a well-defined rises and fall of oscillation except in those cases with large values of 
variance in the photon number such as the left column of section 6 and the left column of section 7. 

Discussion 

In this article, simulated emission and absorption of radiation in a single mode (N-TLMs) for various 
photon distributions has been examined. The original purpose of this paper was to complete (wrap up 
the results from the previous work presented in references L 10 snd 15]. However, as seen below 
and unexpected results was encountered. 

The examples of simulated emission have produced the expected results based on previous work for a 
single TLM for the various photon density distributions and the results for the coherent photon density 
distribution represented by Eq. 17 for various TLM values. On the other hand, the results for simulated 
absorption, especially for the case of small mean photon number compared to the number of TLMs was 
unexpected. (Namely the well-defined rise and fall of oscillation as a function of time. This was 
unexpected since values of photon number approaching or slightly larger than the number of TLMs is 
considerably more chaotic. This may be due to either a coherent effect of the large number of TLMs all 
at equivalent mode positions or due to the photons only interacting with a small number of the TLMs. 


This last conjecture is postulated due to the chaotic effect seen for small photon number when the 
photon variance is large as seen in sections 6 and 7 of Table 4. 

As mentioned above, there is an additional interesting observation concerning simulated absorption 
which occurs when the mean number of photons approaches the number of TLMs. Namely the mean 
number of photons absorbed becomes larger than 34 the number of the initial mean number. This 
behavior was observed for the coherent case [Eq. 17] for all the number TLMs considered above 5 TLMs 
and when fit to a second order equation (the second order component has a very small coefficient) 
showed remarkably small variance from the fit (See Fig. 1 below). 

Note 

There are two included files with this paper. The first is the Wolfram Mathematica notebook (version 
10.0.1) used to provide the calculations and figures presented in this paper. The second is a Computable 
Document Format (CDF) for Interactive Content document. The player can be obtained at Wolfram at no 
expense at https://www.wolfram.com/cdf/. 



Table 2: Photon Density and Si(n,N,yt) 


Coherent Resonant Case 100 Photons 10 TLMs 
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Density Matrix for CoherentState 
for He = nmean=100.0, AVar=100.000 
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Non Resonant Coherent Case 100 Photons 10 TLMs jg=10 or J=25 


Pm CoherentCase 


Density Matrix for CoherentState 

fornc = nmean=100.0, nvar=100.000 


Si(100,10,yt) 




Non Resonant Coherent Case 30 Photons 10 TLMs 0=20 or J=100 
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Table 2: Photon Density and 5i(n,jV,yt) 


Non-Resonant Mixed Coherent Plus Thermal State: 100 Coherent Photons, 30 Thermal Photons and jg=20 or J=100 


Density Matrix Corienent Pluts Thermal State 
torSf=30,0a, 

1130.0. Oyij»7ffl3Q.0a 
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Resonant Squeezed Coherent case 100 coherent photons, squeezing parameter r=l, phase - n and 10 TLMs 


Erenakty Matrix For die Squeezed Coneremt State 
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artel plia*fre3,142 
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Non-resonant Squeezed Thermal State with 100 thermal photons, squeezing parameter r= 1.522, lOTLMs and 0=20 or J=100 


Density Mcitrix for the SqueezedThermal State 
with squeezingparameterr =1.522, Thermal Photon Number= 100.0 
nmean=1056.59, rVar=2.2248 E6 and Critical SqueezingParameter=2.652 


Si(1056.59,10,Kt) 



Non-resonant Squeezed Vacuum case, squeezing parameter r=l, 10 TLMs and jg =20 or J=100 


Density Matrix for the Squeezed Vacuum State 
for the squeezingparamerter=1.000 
nmean=1-381, n^=6.5n06 
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Table 2: Photon Density and 5i(n,jV,yt) 


Resonant Squeezed Fock State with 100 Initial Photons, squeezing parameter r=l and 10 TLMs 


Density Matrix for the SqueezedFock Stcite 
withsqueezingparameterr =1.000, Initial Photon Number=100.0 
^mean=377.6, /7\/ar=66434.9 


Pnn 

0.025 




500 ' 


Non-Resonant Squeezed Fock State with 100 Initial Photons, squeezing parameter r=l, 10 TLMs and jg =20 or J=100 


Density Matrix for the SqueezedFock Staite 
withsqueezingparameterr =1.000, Initial Photon Number=100.0 
^mean=377.6, n\/ar=66434.9 



Si (377.6,10,yt) 
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4 
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200 


300 


400 


500 ' 


Displaced Squeezed Thermal State with 5 thermal photons. 100 equivalent coherent photons, squeezing parameter r=l. phase = n and 10 TLMs 


[>efislty Matrix for the DSTS 

for ri]r =S.0CiD, | ft || ^=1 OO'.d Squeezing Fa ra meter r=1.0(lO 
nv^»S74.'Gg3. phBSea3.142 



Si(120.2,10,yt) 
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Table 2: Photon Density and 5i(n,jV,yt) 


Resonant Squeezed Coherent Plus Thermal State with 100 coherent photons. 5 thermal photons, squeezing parameter r=l, phase = ti and 10 TLMs 


Density Matrix for Squeezed Cofierent Plus Thermal State 
for nr sS.OOO, |P|'=100.0 Squeezing Parameter rsi.ooo 
nnM«-106.4, rrv«-1789.29 



100 150 200 250 300 350 


Si( 106 A 10 ,Kt) 


100 
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300 
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500 ’ 


Non-Resonant Squeezed Coherent Plus Thermal State with 100 coherent photons. 


5 thermal photons, squeezing parameter r=l. phase = n.lO TLMs and 0 =20 or 4=100 


Density Matrix for Squeezed Coherent Plus Thermal State 
for nr =5.000. |P|^=100.0 Squeezing Parameter rsi.ooo 
r»v*-17«9.29 



Si(106.4,10,Kt) 

5 



100 


200 


300 


400 


500 ’ 


Resonant Displaced Number state (equivalent to 100 coherent photons) with initial photon number=25 and 10 TLMs 


Density Matrix for the Displace Number State 
with Initial Photon Number ■25.00 and Number of Coherent PhotonsalOO.O 
nm.-,-125.0. r»v.*-5100.00 



Si(125.0,10,v4) 
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Table 2: Photon Density and Si(n,N,yt) 


Non-Resonant Displaced Number state (equivalent to 100 coherent photons) with initial photon number=25, lOTLMs and 0 =20 or J=100 

Si(125.0,10,/t) 

5 ^ 


Density Matrix for the Displace Number State 
with Initial Photon Number ■25.00 and Number of Coherent PhotonsalOO.O 
«oi**n"125.0. r»v^-5100.00 


Pnn 



Resonant Squeezed Displaced Number state 5 initial photons, displacement equivalent to 100 coherent photons, squeezing parameter r=l. phase =ti andlO TLMs 


Density Matrix Tor Oie Squeezed Dispteced Number State 
for Pnitisli Photoo Number =5 .,Cmjq and Number of cohefeni: Photona=l05,0 
with Squeezing Parameter= 1.005 and Total Plhas«= 5.142 
120.2 a nd ^^*1 =0531.05 


Si(120.2,10,Kt) 



3 

2 
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100 


200 


300 


400 


500 ' 



Non-Resonant Squeezed Displaced Number state 5 initial photons, displacement equivalent to 100 coherent photons, squeezing parameter r=l. phase =ti. 10 TLMs and 0 =20 or 4=100 


Density Matrix Tor Oie Squeezed Displaced Num be r State 
for Iniltisli Photon Number =5,500 and Number of cohereht Phefofia=l05,0 
with Squeezing Parameter= 1.005 and Total Plhase= 3.142 
rim^,n=120.2 and: ^^^,=35311.85 


5,(120.2,10,yt) 

2.8 



405 



Non-Resonant Squeezed Displaced Number state 10 initial photons, displacement equivalent to 100 coherent photons, squeezing parameter r=2. phase =ti. 10 TLMs and 0 =20 or 4=100 


Density Matrix Tor the Squeezed Displaced Number State 
iai Photon Number bIO.OO and Number of CoherenI Photon05.0 
with Squoozing ParamcIcT- 2.005 and Total Phase- 3,142 
nm,.,n=5a5.2 a rrdi =155385. 


Si(386.2,10,yt) 

3.8 

3.6 






































































Table 2: Photon Density and 5i(n,jV,yt) 


Non-Resonant Squeezed Displaced Number state 10 initial photons, displacement equivalent to 100 coherent photons, squeezing parameter r=2. phase =0,10 TLMs and 0 =20 or J=100 


0«nslty Matrix for th« Sque«zed Displaced Number State 
for Initial Pboton Number bIO.OO arKl Number of Coberertt Pbotons-100.0 
with Squeezing Parameter^ 2.000 artd Total Phaser 0 
_ 2 and /v,-41371.5 _ 


Phn 

0.015 




200 400 600 800 1000 1200 



Resonant Squeezed Displaced Number state 10 initial photons, displacement equivalent to 100 coherent photons, squeezing parameter r=l. phase =0 and 10 TLMs 


Density Matrix for the Squeezed Displaced Number State 
for Initial Photon Number bIO.OO and Number of Coherent PhotonsslOO.O 
with Squeezing Parameters 1.000 and Total Phases 0 
_ n,^s130.0 and nv,s10l4.26 _ 


Si(139.0,10,Kt) 

10f 


Pnn 
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0 
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: / 

y 

y 

L. 



Non-Resonant Squeezed Displaced Number state 10 initial photons, displacement equivalent to 100 coherent photons, squeezing parameter r=l. phase =0,10 TLMs and 0 =20 or ^=100 


Density Matrix for the Squeezed Displaced Number State 
for Initial Photon Number sio.oo and Number of Coherent Photonssioo.o 
with Squeezing Parameters 1.000 and Total Phases 0 
_ n,.Has.«130.0 and nv,s10l4.26 _ 


Pnn 

0.035 


0.030 

0.025 
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0.015 
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0.005 



Si(139.0,10,Kt) 
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Table 3: Effect of Increasing Photon number on 5i(n, N, yt) for the coherent photon density 


Pm CoherentCase 
0.07 

0.06 
0.05 
0.04 
0.03 
0.02 
0.01 


[Density Matrix for CoherentState 
for 77c = nmBan =30.00, AVar=30.0000 



80 100 120 140 



Pm CoherentCase 
0.07 

0.06 
0.05 
0.04 
0.03 
0.02 
0.01 


Density Matrix for CoherentState 
for Tic = /7mBan =30.00, AVar=30.0000 


Sif30,10,yt) 
10 




20 40 60 80 100 120 140 


Pm CoherentCase 
0.04 


0.03 


0.02 


0.01 


Density Matrix for CoherentState 

for77c = nmBan=100.0, 77var=100.000 



Si(100,1,yt) 

1.0 
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50 100 150 200 250 300 



200 


400 


600 


800 ' 


Pm CoherentCase 
0.04 


0.03 


Density Matrix for CoherentState 
for7ic = nmBan=100.0, 77var=100.000 


Si(100.0,10,yt) 
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Table 3: Effect of Increasing Photon number on 5i(n, N, yt) for the coherent photon density 


Density Matrix for Coherent State 

for he * nmeafi*200.0, />va,>200.000 


Phn Coherent Case 



Si {200.0,1,yt) 

1.0r 


0.& 


0.4 


0.2 





Density Matrix for Coherent State 

for Tic a nmeaii*200.0, /»va,>200.000 


Phn Coherent Case 





Density Matrix for Coherent State 
for 7»c >nmaan» 1000., hya,-1000.00 


Si(1000.,1,yt) 


p„„ Coherent Case 

0.012 
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A A AAAAJ 

f f f f tiri 


Density Matrix for Coherent State 
for he > /Imaan-IOOO., h^a, >1000.00 


Si (1000.,50,yt) 
50 
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Table 4: Photon Density with a comparison between 54 ( 11 , N, yt) and Si(n, N, yt) 


Density Matrix for Coherent State 



Density Matrix for Coherent State 

for tic — ^mean = 10) 



for Hq = nniean~1^^) 

Calculation Time=0 sec 



Calculation Time=0.01560 sec 


Pnn Coherent Case 



Pnn Coherent Case 




5 TLJ\.l3 all up 
ReBonant Case 

'f5rn?Bn=^Cl, = 

Calculation Ttme=11 se 


5 TLMs all up 
Resonant Case 

lmean=100. nvar=100 

Calculation Tlme=75 sec 




Si(100,5,Kt) 

5 



S4(10,5,Kt) 


50 


100 


200 


400 


5 TLMs all down 



5 TLMs all down 

Resonant Case 



Resonant Case 

nmean = 10, /7var = 10 



/l™an = 10a /Ivar = 100 

Calculation Time=9 sec 



Calculation Tinie=75 sec 


S4(100,5,yf) 
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800 ' 
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Table 4 : Photon Density with a comparison between 84(11, N.yt) and Si(n,N, yt) 


Density Matrix for Coherent State 
for tic = = 

Calculation Time=0 sec 


Density Matrix for Coherent State 
for tic - ^mean=25, n\jar=25 
Calculation T1me=0 sec 


Pnn Coherent Case 

0.12 


Pnn Coherent Case 




100 TLMs all up 



100 TLMs all up 

Resonant Case 



Resonant Case 

^mean = 10t /7\/ar = 10 



/^mean=25, /7var=25 

Calculation Time=3142 sec 



Calculation Time=6249 sec 


Si(10,100,yt) 
lOOr 


Si(25,100,yt) 
100 



54 ( 1 0,100,yt) 
10 


100 TLMs all down 



100 TLMs all down 

Resonant Case 



Resonant Case 

nmean = 10. /7var = 10 



^mean=25, /7var=25 

Calculation Time=139 sec 



Calculation Time=951 sec 


S4(25,100,yt) 



25 
















































Table 4: Photon Density with a comparison between 84 ( 11 , N.yt) and Si(n,N, yt) 


Density Matrix for Coherent State 



Density Matrix for Coherent State 

for tic = nnfiean=1) ^ar=1 



for Tic = nmean = 100, /)var=100 

Calculation Time=0 sec 



Calculation Time=0.01560 sec 


Pnn Coherent Case 



Pnn Coherent Case 




100TLMS all up 
Resonant Case 

^mean = 1> = 1 

Calculation Time=822 sec 


100TLMS all up 
Resonant Case 

^mean“^^^' 

Calculation Time=8504 sec 


Si (1,100,Kt) 
80 


Si(100,100,Kt) 

lool 


3 



S4(1,100,Kt) 


50 


100 


100 


100 TLMs all down 


100 TLMs all down 

Resonant Case 


Resonant Case 

^mean = 1. ^ar = 1 


^mean”^^^’ 

Calculation Time=3 sec 


Calculation Time=10900 sec 


54 ( 1 00,100,Kt) 


150 


200 




Kt 


Kt 














































Table 4 : Photon Density with a comparison between 84(11, N.yt) and Si(n,N, yt) 

Densitv Matrix for Thermal State Densitv Matrix for Thermal St; 


Density Matrix for Thermal State 
for /tj- =nmean=”l> 

Calculation Time=0 sec 


Density Matrix for Thermal State 

for nr =nmean=10) ft\/ar =110 

Calculation Time=0 sec 


Phn Thermal Case 
0.5 


Pnn Thermal Case 



0.04 



6 8 10 12 14 " 

20 40 60 

5 TLMs ail up 

Resonant case 

^fn£arv=1’ 

Calculation Tlnne=3 sec 



lOOTLMsall up 
Resonant Case 

Calculation Time=6069 sec 



Si(10,100,yt) 


20 



50 

100 

150 200 


50 

100 


5 TLMa all down 

Resonant case 

^mean=1r /lvar=2 

Calculation T[me=2 sec 



100 TLMs all down 

Resonant Case 

^mean=10. nvar=110 

Calculation Time=5054 sec 


S4{1,5,yt) 



54 ( 10 , 1 00, yt) 
10 








































































































































































































































Table 4 : Photon Density with a comparison between 84(11, N.yt) and Si(n,N, yt) 


Pnn 


Density Matrix for Coherent Pius Thermai State 



Density Matrix for the Squeezed Coherent State 

for Hr =1, nc=5 



for the squeezing parameter r =1. 

^mean=6) = 



n,=5 

Calculation Time=0.07800 sec 



«mean=6, «var=7 and phase=3.142 



Calculation Time=0.06240 sec 




20 


30 


Si (6,50,/t) 


50 TLMs all up 
Resonant Case 

^mean=6, rt\jax=l 

Calculation Time=549 sec 


5 



50 TLMs all down 



50 TLMs all down 

Resonant Case 



Resonant Case 

^mean=6, nvar^l^ 



^mean=6, /7var=7 

Calculation Tinne=58 sec 



Calculation Time=174 sec 


84 ( 6 , 50 , Kt) 













































Table 4: Photon Density with a comparison 


Density Matrix for the Squeezed Thermal State 



Density Matrix for the Squeezed Vacuum State 

with squeezing parameter r =1.000, Thermal Photon Number=5 



for the squeezing paramerter =1 

^mean=20, /7var=825 and Critical Squeezing Parameter=1.199 



^mean“1> ^ar“® 

Calculation Time=1.076 sec 



Calculation Tlme=0.07800 sec 


Pnn 

0.141 


0.12 

0.10 

0.08 

0.06 

0.04 

0.02 


Pnn 

0.6 

0.5 




50 TLMs all up 



50 TLMs all up 

Resonant Case 



Resonant Case 

^mean=20, /7var=825 



^mean = 1. ^Var^® 

Calculation Time=1562 sec 



Calculation Time=633 sec 


Si(1,50,yt) 





50 TLMs all down 



50 TLMs all down 

Resonant Case 



Resonant Case 

^mean-20, /7var=825 



^mean = 1< ^ar=® 

Calculation Time=1226 sec 



Calculation Time=245 sec 


S4(1,50,yt) 



54(20,50,Kt) 
15 





































































Table 4: Photon Density with a comparison between 84 ( 11 , N.yt) and Si(n,N, yt) 


Density Matrix for the Squeezed Fock State 



Density Matrix for the DSTS 

with squeezing parameter r =1, Initial Photon Number=5 



for Hj =1.000, 1 p 1 ^=5.000 Squeezing Parameter r=1.000 

^mean=20, /7var=203 



^mean = 'IO. '^/ar=63, phase=3.142 

Calculation Time=0.6084 sec 



Calculation Time=0.03120 sec 


Pnn 





40 60 ^ 

10 20 30 

50 TLMs all up 

Resonant Case 

^mean=20, /7var=203 

Calculation Time=1260 sec 



50 TLMs all up 

Resonant Case 

^mean = 10, nvar=63 

Calculation Tinne=726 sec 



Si(10,50,Kt) 


20 


10 



100 


200 


300 


400 


S 4 (20,50, 
20 


15 


50 TLMs all down 



50 TLMs all down 

Reson ant Case 



Resonant Case 

^iaiL‘iin“20, /Jyj,=20'3 



^mean = 10. War =63 

Calculation Tirfie=894 sec 



Calculation Tinne=298 sec 





















































Table 4 : Photon Density with a comparison between 84(11, N.yt) and Si(n,N, yt) 



Density Matrix for Squeezed Coherent Plus Thermal State 



Density Matrix for the Displace Number State 

for Ht =1.000, 1 p 1 ^=5.000, Squeezing Parameter r=1.000 



with Initial Photon Number =1 and Number of Coherent Photons=5 

/’mean=7, /7var=58 and Phase= 3.142 



^mean=6, /7var=15 

Calculation T1me=0.1872 sec 



Calculation Time=0 sec 


Pnn 



50 


100 


150 


50 TLMs all up 



50 TLMs all up 

Resonant Case 



Resonant Case 

^mean=7, n\/ar=^58 



^mean=6. nvar = 15 

Calculation Time=1993 sec 



Calculation Time=263 sec 



Si(7,50,Kt) 


Si(6,50,Kt) 
40 



100 


200 


300 


400 ' 


50 TLMs ail down 



50 TLMs all down 

Resonant Case 



Resonant Case 

^mean=7, /?Var=58 



^mean-6. 

Calculation Time=1631 sec 



Calculation Time=19 sec 


S4(6,50,Kt) 



S4(7,50,Kt) 














































Table 4: Photon Density with a comparison between 84 ( 11 , N.yt) and Si(n,N, yt) 


Density Matrix for the Squeezed Displaced Number State 
for Initial Photon Number =1 and Number of Coherent Photons=5 
with Squeezing Parameter= 1 and Total Phase= 0 

^mean“^0 ^ar“21 

Calculation Tlme=0.01560 sec 


Pnn 

0.20 



20 


40 


60 


100 


Si(10,50,yt) 


50 TLMs all up 
Resonant Case 

^mean”^^' ^ar“21 

Calculation Time=1237 sec 



50 TLMs all down 
Resonant Case 

^mean = 10> /7var=21 

Calculation Time=859 sec 
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